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Solution for a stationary spherically symmetric accretion of the relativistic perfect fluid with an equation of 
state p(p) onto the Schwarzschild black hole is presented. This solution is a generalization of Michel solution 
and applicable to the problem of dark energy accretion. It is shown that accretion of phantom energy is 
accompanied with the gradual decrease of the black hole mass. Masses of all black holes tend to zero in the 
CD phantom energy universe approaching to the Big Rip. 

o : 

_i 1 Introduction 

£ : 

1—5 ■ Our Universe is seems to undergo a period of accelerated expansion and it is assumed that a 
. considerable part of the total density consists of dark energy component with negative pressure pQ. 
There are several candidates for the dark energy: cosmological constant (A) or dynamical component 
7—1 ■ such as quintessence |2J and A;-essence [3]. In connection with the solving of the problem of fine- 
. tuning the models of dynamical dark energy component are seem to be more realistic as they admit 
t— I | to construct "tracker" [I] or "attractor" f.3 ; solutions. 

One of the peculiar feature of the cosmological dark energy is a possibility of the Big Rip [3]: 
■ the infinite expansion of the universe during a finite time. The Big Rip scenario is realized if a dark 
! energy is in the form of the phantom energy with p + p < 0. In this case the scenario of Big Rip 
is possible when cosmological phantom energy density grows at large times and disrupts finally all 
^ | bounded objects up to subnuclear scale. Note, however that the only condition p + p < is not 
i enough for the realization of Big Rip In [7] the authors analyzed the supernova data in the model 
^jrj| independent manner and showed that the presence of the phantom energy with —1.2 < w < — 1 
is preferable in the present moment of time. The analogy between phantom and QFT in curved 
"-j . space-time has been developed in (HJ. The entropy of the universe with phantom energy is discussed 
rS ' in El. 

Usually the evolution of quintessence or /c-essence are considered in a view of cosmological 
problems. However in the presence of compact objects such as black holes the evolution of dark 
energy should be sufficiently different from that in the cosmological consideration. Indeed, what 
would be the fate of black holes in the universe filled with the phantom energy and coming to Big 
Rip? Recently we showed that all black holes gradually decrease their masses and very near the 
Big Rip they finally disappear |lUj . In the present work we study in details the stationary accretion 
of dynamical dark energy into the black hole. As a model of DE we take the perfect fluid with 
negative pressure. The studying of accretion of perfect fluid on the compact objects originated from 
Bondi The relativistic generalization of the perfect fluid accretion were made by Michel |12j . 
Below we find the solution for a stationary accretion of the relativistic perfect fluid with an arbitrary 
equation of state p(p) onto the Schwarzschild black hole. Using this solution we show that the black 
hole mass diminishes by accretion of the phantom energy. Masses of all black holes gradually tend 
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to zero in the phantom energy universe approaching to the Big Rip. The diminishing of a black hole 
mass is caused by the violation of the energy domination condition p + p > which is a principal 
assumption of the classical black hole 'non-diminishing' theorems ^3]. The another consequence of 
the existence of a phantom energy is a possibility of traversable wormholes |14j . In |15 | \T7 \ \1H\ 
authors studied the accretion of scalar quintessence field into the black hole, using the specific 
quintessence potentials V(4>) for the obtaining of the analytical solution for the black hole mass 
evolution. We use essentially different approach for the description of DE accretion into black hole, 
namely, we model the DE by the perfect fluid with the negative pressure. 



2 General equations 

Let us consider the spherical accretion of dark energy onto black hole. We assume that the density of 
the dark energy is sufficiently low so that the metric can be described by Schwarzschild metric. We 
model the dark energy by a perfect fluid with energy-momentum tensor: T^ u = (p + p)u^u u — pg^u, 
where p is the density and p is the pressure of the dark energy and is the four-velocity = 
dx^/ds. The integration of the time component of the energy-momentum conservation law T^\ v = 
gives the first integral of motion 

/ 2 \ 1/2 

(p + p) ( 1 - - + u 2 J x 2 u = C 1 , (1) 

where x = r/M, u = dr/ds and C\ is a constant determined below. 

Given the equation of state p = p(p), one can introduce the function n by the relation: 

dp dn ^ 



p + p n 

The function n plays the role of concentration of the particles, though one can use n for the media 
without introducing any particles. In this case n is the auxiliary function. For general equation of 
state p = p(p), from (J2J) we obtain the following solution for n: 

n (p) ( f d P' ! 

" exp \J7+FF) ] - (3) 



From the conservation of energy- momentum along the velocity UpT^ u - v = 0. using © we obtain 
the another first integral: 

n ( n) _ 

-ux 2 = -A, (4) 



n r 



where rioo (the concentration of the dark energy at the infinity) was introduced for convenience. In 
the case of inflow u = (dr/ds) < and the constant A > 0. From Q and (@J one can easily obtain: 

P + P 2 . 2V /2 n /r\ 

1 Yu = C 2 , (5) 



n \ x 
where 

n Poo +P(Poo) 

n(poo) 

We will see below that the constant A which determines the flux is fixed for fluids with dp/ dp > 0. 
This can be done through finding of the critical point. Following Michel |12j we obtain the para- 
meters of critical point: 

u 2 = — V 2 = ^ (7) 
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where 

y2= n d( P + P )_ 1 
p + p dn 

From this by using © it follows that V 2 = c 2 (p), where c 2 = dp/ dp is the squared effective speed 
of sound in the media. Combining the Eqs. ©, ©, and © we find the following relation: 

p* +p(p*) r ,32/ M 1/2 Poo+_P(poo) , Q v 

which gives the p* for arbitrary equation of state p = p(p). Given p* one can find n(p*) using © 
and values x*, it*, using (j7J) and ©. Then substituting the calculated values in © one can find 
the constant A. Note that there is no critical point outside the black hole horizon (x* > 1) for 
c 2 s < or c 2 s > 1. This means that for unstable perfect fluid with c 2 < or c 2 > 1 a dark 
energy flux onto the black hole depends on the initial conditions. This result has a simple physical 
interpretation: the accreting fluid has the critical point if its velocity increases from subsonic to 
trans-sonic values. In a fluid with a negative c 2 or with c 2 > 1 the fluid velocity never crosses such 
a point. It should be stressed, however, that fluids with c 2 < are hydro dynamically unstable (see 
discussion in [201 G^)- The Eq. © together with © and ((3J) describe the requested accretion flow 
onto the black hole. These equations are valid for perfect fluid with an arbitrary equation of state 
p = p(p), in particular, for a gas with zero-rest-mass particles (thermal radiation) and for a gas 
with nonzero-rest-mass particles. For a nonzero-rest-mass gas the couple of equations © and © 
is reduced to similar ones found by Michel ^2]- One would note that the set of equations ©, © 
and © are also correct in the case of dark energy and phantom energy p + p < 0. In this case 
concentration n{p) is positive for any p and constant C2 in © is negative. 

The black hole mass changes at a rate M = — 47rr 2 T r due to the fluid accretion. With the help 
of JIJ) and © this can be expressed as 

M = 47rAM 2 [ Poo +p{p OD )}. (10) 

For the phantom energy the relation IjlflJI leads to the diminishing of the black hole mass. That 
means that in the universe filled with phantom energy the black holes should melt away. This result 
is general, it does not depend on the equation of state p = p(p), the only condition p + p < is 
important. 



3 The analytical models 

Let us consider the model of dark energy with linear dependence of pressure from the density: 

p = a(p-p ), (11) 

which include, among others, the ultra-relativistic gas (p = p/3) and simplest models of dark energy 
(po = and a < 0). Introduced value a is connected with usual equation of state w = p/p by the 
relation if = a(p — po)/p- For a < there is no critical point for the flux of the fluid into the black 
hole. In the case of a > 0, using (jJJ) and © we find the parameters for critical point in model (|llj) : 

1 + 3a 2 a 
X * = ^a-' U * = TT3a ' (12) 

It should be noted that in the linear model (jlljl the parameters of critical point (|12|) determined 
only by dp/ dp = a and do not depend on the parameter po, which determines what physical fluid 
is considered: relativistic gas, dark energy or phantom energy. Note also that for a > 1 (that 
corresponds to the non-physical situation of superluminal speed of sound) there is no critical point 
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outside the black hole. Let us calculate the constant A which determines the flux of the fluid into 
the black hole. From Eq. © we find: 



n 



Peff 

PeS,cc 



(13) 



where we defined the effective density p e $ = p + p = —poa + (1 + a)p. Using © we obtain: 

a/(l+a) 



PeS,oc 



Jl + 3a) 



1/2 



(14) 



where p e s t is the value of effective density at the critical point and /? e ff,oo is the effective density at 
the infinity. Substituting ((11)) in ((T3*|) and then using @ we find for linear model: 



(1 + 3a)( 1+3 ")/ 2a 

4a 3/2 



(15) 



It is easily seen that A > 4 for < a < 1. For a = 1 that corresponds to c s = 1 the we have A = 4. 
From this we may conclude that for typical sound speeds the constant A has value around unity. 
For some particular choices of parameter a the values p(x) and u(x) can be calculated analytically. 
For example, for a = 1/3 the fluid density is given by: 



Po , ( Po 



z + 



i 



3(1 - 2X- 1 ) 



(16) 



where 



2w _ §_ 

3 3 ; ' 

3 J ' 



x > 3, 



/5 = cos 



and 



1 



2 (a/3) 3 / 2 
2 

2^3 

si 



108 



3(1 -ir 27(1 -r i 1 

The density distribution for another physically interesting case a = 1 is given by: 

2\ / 4 



4 



Po 
4 



1 + 



1 



(17) 



The corresponding radial fluid velocity u = u(x) can be calculated by substituting of (|T6|) or (|T7|) 
into (pQ). For po = the solutions Q16JI and (|17|) describe correspondingly a thermal radiation and 
a fluid with ultra-hard equation of state. In the case of p^ < apo/(l + a) the solutions (|16j) and 
l|17|l describe the phantom energy falling onto the black hole. For example, a phantom energy 
flow with parameters a = 1 and po = ^Poo results in a black hole mass diminishing with the rate 
M = -8ir(2M) 2 Poo . 



4 Black holes in the universe with Big Rip 

Now we turn to the problem of the black hole evolution in the universe with the Big Rip when a 
scale factor a(t) diverges at finite time For simplicity we will take into account only dark energy 
and will disregard all others forms of energy. The Big Rip solution is realized for in the linear 
model (|TT|) for p + p < and a < — 1. From the Friedman equations for the linear equation of state 
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model one can obtain: \p + p\ oc a 3 ( 1+a ). Taking for simplicity po = we find the evolution of the 
density of a phantom energy in the universe: 

where 

_! 3(1 + a) /8vr \ 1 2 



Poo,ij (19) 

and Poq^ is the initial density of the cosmological phantom energy and the initial moment of time 
is chosen so that the 'doomsday' comes at time r. From (|18|) and ()19|) it is easy to see that the Big 
Rip solutionis realized for a = dp/ dp < —1. In general, the satisfying the condition p + p < is not 
enough for the possibility for Universe to come to Big Rip. From (|10(l using (|18f) we find the black 
hole mass evolution in the universe coming to the Big Rip: 

M = M i (l + J^—) , (20) 

V M TT-tJ 1 > 

where 

M = (3/2)A- 1 |l + a|, (21) 

and Mi is the initial mass of the black hole. For a = —2 and typical value of A = 4 (corresponding 
to itH = — 1) we have Mq = 3/8. In the limit t — > r (i.e. near the Big Rip) the dependence of black 
hole mass on t becomes linear, M ~ Mq (t — t). While t approaches to r the rate of black hole 
mass decrease does not depend on both an initial black hole mass and the density of the phantom 
energy: M ~ —Mq. In other words masses of all black holes in the universe tend to be equal near 
the Big Rip. This means that the phantom energy accretion prevails over the Hawking radiation 
until the mass of black hole is the Planck mass. However, formally all black holes in the universe 
evaporate completely at Planck time before the Big Rip due to Hawking radiation. 



5 Scalar field accretion 

In remaining let us confront our results with the calculations of (not phantom) scalar field accretion 
onto the black hole ^3 El QZ| QH| • The dark energy is usually modelled by a scalar field (ft with 
potential V(<p). The perfect fluid approach is more rough because for given 'perfect fluid variables' 
p and p one can not restore the 'scalar field variables' (ft and V</>. In spite of the pointed difference 
between a scalar field and a perfect fluid we show below that our results are in a very good agreement 
with the corresponding calculations of a scalar field accretion onto the black hole. 

The Lagrangian of a scalar field is L = K — V, where K is a kinetic term of a scalar field (ft 
and V is a potential. For the standard choice of a kinetic term K = <ft ;ll (ft'^/2 the energy flux is 
To r = 4> t t<ft,r- Jacobson |15j found the scalar field solution in Schwarzschild metric for the case of 
zero potential V = 0: (ft = (ftoo\t + 2Mln(l — 2M/r)], where (ftoo is the value of the scalar field at 
the infinity. In it was shown that this solution remains valid also for a rather general form 
of runaway potential V(cft). For this solution we have T r = -(2M) 2 ^/r 2 and correspondingly 
M = 4vr(2M) 2 ^. 

The energy-momentum tensor constructed from Jacobson solution completely coincides with one 
for perfect fluid in the case of ultra-hard equation of state p = p under the replacement — > 4^/2, 
Poo — > (ftoa/'Z- It is n °t surprising because the theory of a scalar field with zero potential V(<ft) is 
identical to perfect fluid consideration [22]. In a view of this coincidence it is easily to see the 
agreement of our result ()10|) for M in the case of p = p and the corresponding result of |151 117j . 

To describe the phantom energy the Lagrangian of a scalar field must have a negative kinetic 
term J5 a , for example, K = — <ft- fi (ft' tJ '/2 (for the more general case of the negative kinetic term see |19j). 
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In this case the phantom energy flux onto black hole has the opposite sign, Tq t = —<j)t<f> r -> where 4> 
is the solution of the same Klein-Gordon equation as in the case of standard scalar field, however 
with the replacement V — > — V. For zero potential this solution coincides with that obtained by 
Jacobson 15 for a scalar field with the positive kinetic term. Lagrangian with negative kinetic term 
and V(4>) = does not describe, however, the phantom energy. At the same time, the solution for 
scalar field with potential V(4>) = is the same as with a positive constant potential Vo = const, 
which can be chosen so that p = — <p 2 /2 + Vq > 0. In this case the scalar field represents the required 
accreting phantom energy p > and p < —p and provides the decrease of black hole mass with the 
rate M = -4^(2M) 2 0^. 

The simple example of phantom cosmology (without a Big Rip) is realized for a scalar field 
with the potential V = m?(j) 2 /2, where m ~ 10~ 33 eV After short transition phase this 

cosmological model tends to the asymptotic state with H ~ m</>/3 1//2 and <j) ~ 2m/3 1 / 2 . In the 
Klein-Gordon equation the m 2 term (with the mentioned replacement V — > — V) is comparable 
to other terms only at the cosmological horizon distance. This means that the Jacobson solu- 
tion is valid for this case also. Calculating the corresponding energy flux one can easily obtain 
M = -4vr(2M) 2 2 X) = -64M 2 m 2 /3. For M = M© and m = 1CT 33 eV the effective time of black 
hole mass decrease is r = (3/64)M _1 m -2 ~ 10 32 yr. 
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